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Abstract — -A numerical reconstruction of the brightness 
temperature is examined as a potential way to improve the 
retrieval of salinity from Aquarius measurements closer to land- 
water boundaries. A test case using simulated ocean-land scenes 
suggest promise for the technique. 
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I. Introduction 

The goal of the Aquarius [1] is to map the surface salinity 
field of the global oceans with a spatial resolution of 150 km 
and a monthly RMS accuracy of 0.2 psu [2], Because of the 
accuracy requirements for the salinity retrieval and the large 
footprint of the Aquarius radiometers, the retrieval of salinity is 
restricted to the open ocean several hundred kilometers from 
land. At closer distances, land in the antenna side lobes, biases 
the measurements and causes an underestimation of the salinity 
field. In principle, it is possible to correct for this 
“contamination” by land if the antenna pattern and emission 
from land are well known. However, this is generally not the 
case. In this study, an alternative approach for retrieving 
salinity close to land is developed using a regularized 
reconstruction of the brightness temperature. 

II. AQUARIUS SALINIY RETRIEVAL 

Aquarius radiometers report measurements of the 
microwave emission from Earth every 1 .44 seconds. The basic 
observation geometry is shown in Fig. 1. After making 
corrections for interfering sources of radiation, such as the 
galactic background, the Sun and the atmosphere [3], the 
antenna temperature at a point along the boresight can be 
written in the form [4]: 

T A( r i)= J G(co,b 0 )T B (co) dco (1) 

B FOVi 

where T A = [T ah ,T av ] t , T B = [T h ,T v ] t and the antenna 
pattern matrix pointing in the boresight direction b 0 , is 
G (a),b 0 ). Here H F0F . is the solid angle subtended by the 
visible disk (FOV) corresponding to the ith measurement. 



Figure 2 illustrates the problem with land in the field-of- 
view of the radiometer. The error due to radiation from land is 
plotted (ordinate) as a function of the distance of the antenna 
boresight from the land-water boundary (abscissa). In this 
example, the outer-most Aquarius radiometer beam is used 
(local incidence angle of 45.6 degrees). An error of 0.1 K 
corresponds roughly to the goal of 0.2 psu. Hence, the 
acceptable error is reached at about 450 km from the land- 
water boundary which is about 4.5 footprints (3 dB) for the 
outer beam. 



Figure 2: The effect of land contamination on the observed 
brightness temperature. 
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The goal here is to retrieve salinity closer to shore without 
being critically dependent on an accurate model for land 
emission to correct for land contribution in the side lobes. 

III. BRIGHTNESS TEMP RECONSTRUCTION 

The alternative approach investigated here is to take 
advantage of the oversampling in the along-track dimension 
inherent in the Aquarius data collection. This is done by 
converting (1) to a discrete integral equation and attempting to 
invert it. The integral equation is ill-posed but can be inverted 
by regularization. The advantage of this approach is that it does 
not require an accurate model for emission from land in the 
FOV. 


The first step is to rewrite (1) in terms of the reference grid 
shown in Fig 3. The box is chosen to be about 2.5 times the 3 
dB width and several beam width long. In order to simplify the 
mapping, the brightness temperature is expanded in a Taylor 
series with respect to incidence angle about boresight. 
Keeping first order terms, one obtains: 


T A (Xi) * J G p(x-ri,b 0 )T B (x,b 0 )dx 

G 

+ J G P (x -r^bo^Cid + c 2 ti 2 ]dx 


( 2 ) 


where d is the difference between the incidence angle at point 
x and the incidence angle, d 0 at boresight point r^, and the a 
are the derivatives of the brightness temperature evaluated at 
boresight. Since the angular dependence of the brightness 
temperature is only weakly dependent on the value of Tb itself, 
and approximate value can be used to compute these 
coefficients. 

Assuming the antenna pattern, G P , and the c’s are known, 
the second term is known and the integral equation (1) can be 
transformed to the convolution below: 

T Ao { r i) » J G P (x — ri,b 0 )T Bo (x ) dx (3) 

G 

where T A o is T A minus the second integral in the previous 
expression and T Bo = T B (x,b 0 ) is the unknown brightness 
temperature field at the fixed incidence angle i9 0 . Solutions to 
(3) are sought using a uniform step quadrature to transform (3) 
to the linear algebraic system: 

Tao = KT Bo ( 4 ) 

where K is the matrix corresponding to the quadrature in (3) 

and T Bo is the column vector of unknown brightness 
temperatures at the discrete quadrature grid points of the region 
G shown in Fig.3. 

There are N x grid points along track and N y grid points 

cross track. Then T Bo is a vector of 2 N x N y elements (the 

factor of 2 is for the two polarizations). Assuming there are N m 
measurements of the antenna temperature for each polarization, 


then T Ao is a vector of 2 N m elements. Since N m < N x N y , the 
linear system in (3) is rank deficient. 





Figure 3: Study region with discrete grid 


The brightness temperature T Bo is reconstructed with the 
help of assumptions which provide the additional equations 
required to regularize the system in (3). The following 
assumptions are made about the solution: 


• Smoothness along track, 

• Flatness cross track (no cross track information) 

• The polarizations are dependent (polarization 
consistency) 

• Proximity to some initial estimate T BE 


With these assumptions the regularized solution to the 
system in (4) is given by: 


Tbo = K'T Ao + L- 1 E t ET be 

where 


( 5 ) 


= L~ x K t (6a) 

and 

L = K t K + E t E + C T C 

(6b) 

+ OC x D 2,X^2,X + ttyD\yDiy 

Here, C is a diagonal matrix with entries a c over water to 
provide polarization consistency, D 2x is a 2 nd order finite 
difference matrix for smoothness along the x-direction (along 
track), D l y is a 1 st order finite difference matrix for flatness 
along the y-direction (across track) and E is a diagonal matrix 
with entries a t and a w over land and water, respectively, for 
proximity to the estimate T BE . 
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The reconstructed solution in (4), requires five 
regularization parameters. These are namely, a c , a t , a W9 a x 
and a y . In this study, they are determined manually by trial and 
error. 

IV. RECONSTRUCTION EXAMPLE 

Example cases have been simulated to demonstrate the 
feasibility of the proposed method. In the examples the grid 
width is +10° which contains slightly more than 98% of total 
power in the Aquarius antenna beam; and the grid length is 
about 1131 km corresponding to 101 of the 1.44-sec data 
points (N m = 101). The grid is divided into N x = 101 by 
N y = 11 points for a total of 1111 grid points. The Aquarius 
outer beam is used in the examples. 

The “data” are simulated antenna temperatures using the 
actual Aquarius antenna patterns and models for emission from 
the surface. This permits a way to test the approach because the 
surface brightness temperature is known. The ocean surface 
model uses sea surface temperature from NCEP [5] and salinity 
from HYCOM [6]. The land surface model uses NCEP values 
for temperature and soil moisture with a vegetation canopy 
using parameters from ECOCLIMAP [7]. 

Finally, the estimate T BE is a constant determined from 
these emission models by adding bias (0.5 psu for salinity, 5% 
for soil moisture and 5 K for land surface temperature) and 
random error to the ancillary data and taking the mean of the 
resulting brightness temperatures for land and sea. 

In the example below, we consider a land crossing in the 
vicinity of the Amazon plume which is a dynamic area where 
large amounts of fresh water pour into the sea and significant 
salinity variations occur over relatively short distances. The 
simulated scene is shown in Fig 4 (left) for V-pol. It was 
generated using ancillary data for August 27, 2011. The 
Amazon River delta is just below the southern edge of the map. 
The dashed line is the beam center (boresight) for the Aquarius 
outer beam. 


T v (K) T y (K) - Reconstructed 



Figure 4: Original scene (left) and reconstructed image 

(right). The simulation is for V-pol and the Aquarius outer 
beam. 

Fig 4 (right) shows the reconstructed image for V-pol. The 
reconstruction along the boresight reproduces the scene quite 
well, although, the quality of the reconstruction degrades as 
one moves to the side away from boresight. This is to be 
expected because the image is two dimensional but Aquarius 


only provides data in one dimension (along track in time). The 
reconstruction has been forced to be uniform across track to 
address this problem (the assumption of flatness across track 
listed above). 

However, along track the reconstruction near boresight is 
very good. This can be seen more clearly in Figure 5 which 
shows the scene brightness temperature along boresight (“x”) 
and the reconstructed brightness temperature (red) as a 
function of the distance from the land-water boundary. The 
error is much less than 0.1K right up to the land-ocean 
boundary. This was done starting with an initial constant bias 
in the estimate for land of 5% for soil moisture and 5 K for 
surface temperature. The blue curves indicate the effect of 
adding 0.2 K of noise to the measurements of antenna 
temperature. The 95% confidence interval (blue curves) are 
within 0. 1 K right up to the land-ocean boundary. 



Figure 5: Performance of the reconstructed brightness 

temperature along boresight as a function of distance from the 
land-water boundary. 

V. CONCLUSIONS 

The test with simulated data suggests the approach for 
retrieving salinity closer to shore is promising. The approach 
provides enhanced spatial resolution, but the goal is to retrieve 
salinity closer to shore. The next step is to test the procedure 
with actual Aquarius data. 
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